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We investigate whether the annual modulation observed in the Dama experiment can be due to a 
weakly-interacting massive particle (WIMP) with an axial- vector (spin-dependent; SD) coupling to 
nuclei. We evaluate the SD WIMP-proton cross section under the assumption that such scattering 
accounts for the Dama modulation, and we do the same for a SD WIMP-neutron cross section. We 
show that SD WIMP-proton scattering is ruled out in a model-independent fashion by null searches 
for energetic neutrinos from WIMP annihilation in the Sun, and that SD WIMP-neutron scattering 
is ruled out for WIMP masses > 20 GeV by the null result with the Dama Xe detector. A SD 
WIMP with mass < 20 GeV is still compatible, but only if the SD WIMP-neutron interaction is 
four orders of magnitude greater than the WIMP-proton interaction. 



PACS numbers: 95.35.+d; 14.80.Ly; 95.55.Vj 
I. INTRODUCTION 

Weakly-interacting massive particles (WIMPs) have 
long been sought experimentally as the primary com- 
ponent of the dark halo that enshrouds the Milky Way 
0. Such dark-matter particles could be detected di- 
rectly via observation of the O(10 keV) nuclear recoils 
they would produce in a low-background detector when 
a dark-matter particle scatters elastically from a nucleus 
therein ||[|. They could also be detected indirectly via 
observation of the energetic neutrinos produced by an- 
nihilation of WIMPs that have accumulated in the Sun 
and/or Earth 

The WIMP can scatter elastically from a nucleus either 
via a scalar (spin- independent; SI) interaction, in which 
case the WIMP-nucleus cross section scales with the mass 
of the nucleus (see, e.g., Fig. 26 in Ref. |J), or via an 
axial-vector (spin-dependent; SD) interaction, in which 
case the cross section is nonzero only if the nucleus has 
a non-vanishing spin. In this latter case, WIMP-nucleus 
scattering occurs primarily via interaction of the WIMP 
with the unpaired proton or neutron, and therefore the 
WIMP-nucleus cross section will mostly depend on the 
SD WIMP-proton interaction in odd-Z nuclei and on the 
SD WIMP-neutron interaction in odd-iV nuclei. 

The Dama collaboration || has reported an annual 
modulation in their Nal detector compatible with the 
summer-winter variation in the flux of WIMPs incident 
on the detector due to the motion of the Earth through 
the halo (3). If interpreted as a WIMP with SI interac- 
tions, the effect singles out a region in the plane (WIMP 
mass, WIMP-nucleon cross section), centered at about 
(50 GeV, 7 x 10~ 6 pb), which is largely excluded by the 
null result reported by the Cdms collaboration using a 
detector made of natural germanium ||. Further data 
from more sensitive detectors will allow a firmer state- 
ment on this contradiction. 

An alternative explanation is that the WIMP has a 
dominantly SD coupling with nuclei. In this case, the in- 



compatibility between the Dama and Cdms results dis- 
appears as Na and I are unpaired-proton nuclei, while 
natural germanium contains only ~ 8% of the 73 Ge iso- 
tope which is an unpaired-neutron nucleus. Thus, it is 
conceivable that a WIMP could undergo a SD interaction 
in one detector while remaining invisible to the other. 
Such an explanation has been neglected since the re- 
quired WIMP SD cross sections exceed those expected 
in currently favored supersymmetric models. However, 
given the number of free parameters in the minimal su- 
persymmetric extension of the standard model (MSSM) , 
not to mention the numerous possible alternatives, a 
WIMP with a strong SD coupling to nuclei is certainly 
plausible. 

Here we show that this possibility can be rejected ex- 
perimentally (except for a small unusual corner of param- 
eter space), rather than through theoretical prejudice. 
WIMPs with SD couplings to protons will be captured 
in the Sun and annihilate therein producing energetic 
neutrinos that should be detectable in several existing 
neutrino telescopes. Following a model-independent ap- 
proach we show that current upper limits to the 
neutrino flux rule out the possibility that the Dama sig- 
nal is due to a WIMP with SD couplings to protons. 
We check also that the alternative hypothesis — that the 
modulation signal is due to SD scattering from neutrons 
in Na and I — is, excluded for WIMP masses M x > 20 
GeV by null searches with odd-neutron targets (the best 
limit comes from an experiment with enriched liquid Xe 
performed by the Dama collaboration as well Q). 

Indirect-detection rates have already been calculated 
for a variety of supersymmetric models that fit the Dama 
modulation signal (see Ref. jl(| and references therein). 
Here, however, we carry out a model-independent analy- 
sis that does not rely on any details of the particle-physics 
model. 

In the next Section, we discuss the procedure to derive 
the region in the plane (WIMP mass, WIMP-nucleon 
cross section) compatible with the Dama modulation 
signal both for SI and SD interactions. Section III re- 
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views the model-independent constraints from Ref. |t]] 
to this parameter space in case of SD interactions with 
protons. Section IV shows that the null result with Xe 
strongly constrains the possibility that the Dama mod- 
ulation could be due to a WIMP SD interactions with 
neutrons. Section V concludes. 



II. ANALYSIS OF THE MODULATION SIGNAL 

Dama has analyzed the measured modulation assum- 
ing it is due to WIMPs with SI couplings only, publishing 
the corresponding compatible region of WIMP masses 
and cross sections. The analysis for WIMPs with SD in- 
teractions is completely analogous. We now detail the 
calculation. 

The differential direct-detection rate (per unit detector 
mass) in a detector made of nucleus i is 
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where Q is the energy deposited in the detector, and 
da X i/dQ is the differential cross section for WIMP elas- 
tic scattering with the target nucleus. We assume here 
that WIMPs of mass M x account for a local dark-matter 
density p x and have a local (i.e., in the rest frame of the 
detector) distribution in velocity space / with normaliza- 
tion / d 3 v'f(v') = 1. 

The differential cross section is usually re- written as 
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where Q max = 2 \v'\ 2 /M i: and 
mass between particles j and i. The form- factor suppres- 
sion of the cross section depends on whether the interac- 
tion is SI or SD. The SI form factor Fi(Q) is relatively 
simple to calculate; the estimate for the SD form factor 
Si{Q, dp, a n ) requires more sophisticated nuclear model- 
ing jy] and differs in case of SD WIMP-proton and SD 
WIMP-neutron scattering. 

To compare results obtained with detectors of different 
materials, the cross section for WIMP scattering from 
nucleus i (for SI or SD) is expressed in terms of the 
WIMP-proton (neutron) cross section scaled to zero mo- 
mentum transfer, a X i — (m, xi /m xp )C p ( n )<j p ( n y For SI 
scattering, C^ n) = [Z z f p + (A, - Z A ) /n] 2 //p(„)> where 
Ai and Z\ are the atomic mass number and charge of 
nucleus i, and f p ~ /„ are the SI couplings of WIMPs 
to protons and neutrons, respectively. For SD scatter- 
ing, Cg£) = 4 (A; (n) ) 2 Ji( J, + I)/3, where J t is the total 
angular momentum of nucleus i, and A*^ = (a p (5p) + 
a n{S l n ))/{a p ^Ji). Here, a p and a n are the WIMP SD 
couplings to protons and neutrons and (<5p) and (S 1 ^) 



are the proton and neutron spin expectation values in 
nucleus i; they are obtained from nuclear models ]Tl[ |. 

Finally, to get the differential detection rate for Nal 
one has to take into account the fact that this is not 
a monatomic material and include quenching factors for 
both Na and I. 

Eq. (Q) is time dependent; the dependence is written 
implicitly in the change of variables between the velocity 
of the WIMP in the detector rest frame and in the galac- 
tic frame. The rotation of the Earth around the Sun gives 
rise to the well-known annual-modulation effect 0. To 
a good approximation, the signal event rate in the fcth 
energy bin can be parameterized by separating a com- 
ponent averaged during the year from a time-dependent 
component, 



Sk = S 0ik + 5 mjfc cos [2ir(t - t )/T] 
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where the period T is I year and t$ is about June 2. 
Dama performs a maximum-likelihood analysis and pro- 
vides values of <So,fc and S m ,k corresponding to their fa- 
vored region in 4 energy bins, with electron equivalent 
(measured) energies Q ee between 2 and 6 keV ||. 

Since we cannot perform an analysis on the raw data, 
we consider an approximate method to decide whether a 
WIMP candidate is compatible with the Dama modula- 
tion or not. We follow Ref. jl2) and define the statistical 
variable, 
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where S^ k p is the experimental result and Sq\ is the the- 
oretical prediction for the WIMP candidate being inves- 
tigated. 

Before considering the SD case of interest to us here, 
we first carry out our analysis with a SI WIMP-nucleon 
coupling to test whether our technique reproduces the 
results of the full Dama analysis that has been published 
for this case. To do so, we evaluate Sq\ and S^ k as a 
function of WIMP mass and cr^ 1 (we assume, as Dama 
did, that f p = f n ', we take also their choice of form factors 
and astrophysical parameters) . The 3er region singled out 
by Dama with their maximum-likelihood method in the 
(M x , dp 1 ) plane is reproduced fairly well if we require 
that k < 60. 

We now turn to our analysis of a WIMP with SD in- 
teractions. We take the form factor for SD scattering 
from Na and I from Ref. |r| and again assume that the 
WIMP candidate is compatible with the Dama modu- 
lation if k < 60. We consider separately the case of 
coupling with protons only (a p ^ 0, a n = 0) and with 
neutrons only (a p = 0, a n ^ 0); the results are indicated 
by the shaded regions, respectively, in Fig. ^ and in Fig. 
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FIG. 1. The shaded region shows the parameter space (in 
WIMP mass versus SD WIMP-proton cross section) implied 
by the Dama annual modulation for a WIMP with exclu- 
sively SD interactions with protons and no interaction with 
neutrons. The solid curve indicates the upper bound to the 
SD WIMP-proton cross section from null searches for neu- 
trino-induced upward muons from the Sun; thus the cross 
hatched region is excluded. 



III. A WIMP-PROTON INTERACTION? 

We examine first the case of WIMPs with couplings 
with protons. Such WIMPs scatter efficiently from pro- 
tons in the Sun, losing enough energy to become grav- 
itationally bound to the Sun. This process creates an 
enhancement in the density of WIMPs at the center of 
the Sun. The WIMPs can then annihilate in the Sun, 
and among the decay products of the annihilation prod- 
ucts will be energetic neutrinos (energies of order half the 
WIMP mass) that can escape the Sun. Such neutrinos 
can produce upward muons via a charged-current inter- 
action in the material below a neutrino detector, such 
as that at 1MB, Super-Kamiokande, Baksan, MACRO, 
and/or AMANDA. 

The calculation of the flux of such neutrino-induced 
upward muons from the Sun is straightforward. For the 
masses we are considering, the flux can be written 

T = 0.016m-V _1 (M x /GeV)(a^ /10" 40 cm 2 ) 

( Px /0.3GeVcm- 3 )S(m x /m H )am x ), (5) 

Here, S(x) is given by Eq. (9.21) in Ref. 0. The func- 



*Note in Eq. (9.55) in Ref . Q, the factor tanh(t0/ro) should 
be squared and there should be a factor of £( m x) on the right- 
hand side. 



tion £(to x ) is given in Fig. 33 in Ref. 0; it quantifies 
the number of neutrino-induced muons expected per an- 
nihilation event. This depends in detail on the annihila- 
tion products; e.g., bb or cc quarks, and/or t + t~ lepton 
pairs. Different WIMP candidates in the mass ranges we 
are considering (m x < 80 GeV) will produce different 
branching fractions to these annihilation products, and 
this will result in some allowable range for £(m x ). For the 
WIMP masses we are considering, 0.03 < £( TO x) 0.13; 
we use £ = 0.03. The bound we compute here could be 
evaded if the WIMP annihilated to uu, dd, ss, e + e~, 
and/or pairs but not cc, bb, nor t + t~ pairs. How- 

ever, we know of no models in which (nor any reason why) 
this would occur. Eq. (refeqmindirectrate) assumes that 
capture of WIMPs in the Sun from the Galactic halo is 
in equilibrium with their depletion by annihilation. As 
argued in Ref. [|lj, capture and annihilation will be in 
equilibrium in just about any model in which the signal 
is anywhere close to being detectable. Since we are going 
to place upper limits to a WIMP-proton scattering cross 
section based on current bounds, we may safely assume 
capture-annihilation equilibration in our analysis. 

The upper limit to the flux of neutrino-induced muons 
from the Sun (roughly < 10~ 2 m~ 2 yr _1 ) from Super- 
Kamiokande [H, leads to the upper limit to the cross 
section for WIMP-proton SD scattering shown in Fig. 
[I] (limits from Baksan, MACRO and AMANDA lead to 
very similar results). The discrepancy of a couple of or- 
ders of magnitude between the region disallowed by neu- 
trino telescopes and the region implied by Dama is much 
too large to be explained by the uncertainties and ap- 
proximations inherent in our analysis, which might con- 
ceivably change our Dama parameter space and/or our 
Super-Kamiokande limit by no more than factors of a 
few. 



IV. A WIMP-NEUTRON INTERACTION? 

In the simplest odd-group model of the Na and I nu- 
clei, the spin with which a SD WIMP interacts is carried 
exclusively by the unpaired proton. In such a model, the 
WIMP can undergo SD scattering with the nucleus only 
through the WIMP-proton interaction, and not through 
any WIMP-neutron interaction. However, in more so- 
phisticated nuclear models, some small fraction of the 
spin is carried by neutrons as well. It is thus conceivable 
that the Dama modulation could have been caused by a 
SD WIMP-neutron interaction. The point of this Section 
will be to show that this possibility is disfavored by the 
null results of WIMP searches with odd-neutron targets. 

In particular, the best current limit has been obtained 
from an experiment with enriched liquid Xe (99.5% 
129 Xe) performed as well by the Dama collaboration (we 
checked that a weaker limit is given by the Cdms mea- 
surement with natural Ge, with only ~ 8% of the 73 Ge 
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FIG. 2. The shaded region shows the parameter space (in 
WIMP mass versus SD WIMP-neutron cross section) implied 
by the Dama annual modulation under the assumption that 
the signal is due exclusively to SD WIMP-neutron scatter- 
ing. The solid curve indicates the upper bound to the SD 
WIMP-neutron cross section from the Dama Xe detector. 
Again the cross hatched region is excluded. 

isotope). We plot in Fig. ^ the limit derived in Ref. ||] 
using again form factors from Ref. p3| , but assuming 
a n = — 0.85 a p rather than a p = 0, o„ ^ 0. The limit in 
the latter case should be even slightly lower than the one 
displayed. The limit from Xe excludes all WIMP masses 
larger than about 20 GeV, leaving just a tiny corner of 
the Dama favored region. 



V. CONCLUSION 

The annual modulation observed in the Dama low- 
energy bins has been attributed to a WIMP with SI in- 
teractions with nuclei, a possibility in contradiction with 
null searches by Cdms. However, this is a priori not the 
only interpretation — the modulation could also be due 
most generally to a WIMP with SD interactions. 

We have shown here that the region of parameter space 
implied by the Dama signal for a WIMP that under- 
goes SD scattering from protons is excluded in a model- 
independent way by null searches for a WIMP-induced 
energetic-neutrino flux from the Sun. 

The alternative hypothesis — that the modulation is 
due to WIMP scattering from neutrons — is excluded for 
M x > 20 GeV by data taken with a Xe detector. Al- 
though experimentally viable, the M x < 20 GeV solution 
requires new physics beyond the minimal supersymmet- 
ric standard model |l5). More importantly, it requires 
that the WIMP-proton and WIMP-neutron interactions 



differ by more than four orders of magnitude. Such a 
large difference would be quite unusual. 

The SD WIMP parameters compatible with the Dama 
modulation derived here may be slightly changed by in- 
cluding particle and nuclear uncertainties or uncertain- 
ties in the local WIMP velocity distribution , as well 
as by performing a more accurate analysis on the raw 
data. However all these effects cannot account for the 
order-of-magnitude (or more) discrepancies pointed out 
here. 

We thank P. Belli and J. Edsjo for useful discussions. 
This work was supported in part by NSF AST-0096023, 
NASA NAG5-8506, and DoE DE-FG03-92-ER40701 and 
DE-FG03-88-ER40397. 



[1 
[2 
[3. 

K 

[5 
[6] 
[7] 
[8] 
[9] 
[10 
[11 
[12 
[13 
[14 

[is; 

[16 



G. Jungman, M. Kamionkowski, and K. Griest, Phys. 
Rep. 267, 195 (1996). 

M. W. Goodman and E. Witten, Phys. Rev. D 31, 3059 
(1986); I. Wasserman, Phys. Rev. D 33, 2071 (1986). 
A. Drukier, K. Freese, and D. Spergel, Phys. Rev. D 33, 
3495 (1986); K. Freese, J. Frieman, and A. Gould, Phys. 
Rev. D 37, 3388 (1988); K. Griest, Phys. Rev. D 37, 
2703 (1988). 

J. Silk, K. A. Olive, and M. Srednicki, Phys. Rev. Lett. 
55, 257 (1985). 

R. Bernabei et al. (Dama Collaboration), Phys. Lett. B 
480, 23 (2000). 

R. Abusaidi et al. (Cdms Collaboration), Phys. Rev. 
Lett. 842000 5699. 

M. Kamionkowski, K. Griest, G. Jungman, and B. 

Sadoulet, Phys. Rev. Lett. 75, 5174 (1995). 

M. Kamionkowski and K. Freese, Phys. Rev. D 55, 1771 

(1997). 

R. Bernabei et al. (Dama Collaboration), Phys. Lett. B 
436, 379 (1998). 

A. Bottino, F. Donato, N. Fornengo, and S. Scopel, Phys. 
Rev. D 62, 056006 (2000). 

J. Engel, S. Pittel, and P. Vogel, Int. J. Mod. Phys. E 1, 
1 (1992). 

M. Brhlik and L. Roszkowski, Phys. Lett. B 464, 303 
(1999). 

M. T. Ressell and D. J. Dean, Phys. Rev. C 56, 535 
(1997). 

A. Okada et al. (Super-Kamiokande collaboration), 
[astro- ph/0007003]. 

G. Abbiendi et al. (OPAL Collaboration), Eur. Phys. J. 
C14, 187 (2000). 

M. Kamionkowski and A. Kinkhabwala, Phys. Re v. D 
57, 3256 (1 998); P. Ullio and M. Kamionkowski, |iei> 
ph/0006183| ; F. Donato, N. Fornengo, and S. Scopel, 



Astropart. Phys. 9, 247 (1998); P. Belli et al., Phys. Rev. 
D 61, 023512 (2000); N. Evans, C. M. Carollo, and P. 
T. de Zeeuw, [astro-ph/0008156]; A. M. Green, [astro- 
ph/0008318]. 



4 



